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ABSTRACT 



Effects of boundary irregularities on acoustical mode 
propagation in ducts, with pressure release walls are 
theoretically studied. Several particular wall irregularities 
are explored: Sinusoidal, 6-function and spatially decaying 

wall perturbations. Strong effects (resonances) appear in 
the case of the sinusoidal wall. In the other two cases, 
the resonances are not present, but there are still traveling 
disturbances whose effects can be important at great 
distances from the wall irregularity. 
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I. 



INTRODUCTION 



The theory of acoustic wave propagation in bounded 
media (waveguides and waveducts) has been treated with a 
certain intensity in the last twenty-five years. It is 
possible to classify the problems arising in these studies 
in two categories. One category involves the effects on 
acoustic propagation of small inhomogeneities and the other, 
the effects on the acoustic field of irregularities on the 
boundaries of the waveguide or waveduct. 

These two factors are the most important concerning the 
acoustic w ave propagation and M.A. Isakovich [4] in 1956 
said that "such problems, besides their theoretical interest 
are also important from a practical point of view; thus, for 
example, in super-long-range waveduct propagation of sound 
in a laminarly-inhomogeneous medium the effect of even small 
inhomogeneities or corrugation of the boundaries can produce 
a very great effect upon the range of propagation" . 

M.A. Isakovich studied propagation in a waveduct with small 
inhomogeneities and in a waveduct with one corrugated wall, 
using the method of small perturbations. This was a natural 
extension of the problem of scattering of sound from one 
single rough wall. 

I.C. Samuels [5], in 1958, studied the problem of 
propagation of acoustic and electromagnetic waves in ducts 
with two corrugated rigid walls. He found, in the case of 
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harmonic roughness , "that waves are not scattered upstream 
if the frequency of the transmitted signal exceeds certain 
critical values related to the natural modes of the waveguide 
and the wavelengths of the wall roughness." 

C.S. Clay [6]/ in 1963 presented a paper where the 
effect of a slightly irregular boundary on waveguide 
propagation was studied. The results were checked by an 
experiment in shallow water and he found that the effects of 
small roughness on the surface are to increase the attenua- 
tion as a function of range and decrease the coherence of 
acoustical propagation. 

In more recent years (1971) , R.H. Ebert [7] calculated 
the fluctuation of sound due to influence of a standing 
gravity wave in a waveguide. The results showed a resonance 
if the wavelength of the surface wave was one-half of the 
acoustic wavelength. In an experiment realized in the Naval 
Postgraduate School he verified the existence of that 
resonance . 

R.F. Salant [8], in 1972, stated that the research and 
study of the theory of sound propagation in waveguides could 
be and had been applied to several practical problems and 
gave as examples: transmission of noise through ventilating 

ducts, jet-engine inlets and other ducting systems. He 
considers propagation in a waveguide having two sinusoidal 
walls having a difference of phase, as opposed by earlier 
work by J.C. Samuels [5], whose waveguide had the sinusoidal 
walls in phase. 
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Salant concluded that the phase difference between the 
walls, the wall wave number and the mode number in the wave- 
guide affect strongly the disturbance generated by the walls. 

He found that the disturbances when odd modes are propa- 
gating in the waveguide are extremely different from the 
ones formed when even modes are propagating. 

A. H. Nayfek [9] in 1974 treated again the acoustic 
propagation in ducts with sinusoidal walls using the method 
of multiple scales, in contrast with previous works where 
the small perturbations method was used. For travelling 
waves, resonance occurs wherever the wall wave number is 
equal to the difference of the wave number of any two duct 
acoustic modes. 

N. E. Davis and M. W. Kenyon [1], in December 1974 
studied the effects of travelling sinusoidal surface waves 
in a waveguide, verifying that resonances occur under certain 
conditions relating wavelengths of surface and acoustic waves. 
A waveguide with pressure release walls was built and the 
presence of such resonances verified. 

In this work, a simple theory is presented for solutions 
of the acoustic velocity potential in irregular ducts with 
pressure release walls using the Fourier transform of the 
functions describing the irregularity of the wall. 
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II. PROBLEM DESCRIPTION 



A. GENERAL FORMALISM 

Consider the existence of a waveduct, with transverse 
dimension & . Let the walls be pressure release and smooth, 
except for the upper wall that will be irregular. This 
waveduct will be excited by a monofrequency sound source, 
and is assumed to be of infinite length in the x direction. 

For the case of perfectly smooth pressure release walls, 
a solution for the wave equation (velocity potential) is, 
taking for simplicity an amplitude of one. 



This represents a wave travelling in the positive x direction. 

In this equation, k and k must obey the relations: 

X z 



Assume that the upper wall of the waveduct is irregular 
and located at 



0 = sin(k z) cos (cot - k x) . 
Z X 



( 2 . 1 ) 



k = [(-) 2 - k 2 ] 1/2 

X c z 



( 2 . 2 ) 




, . . . 



z 



z = £ z + e f(x) 



( 2 . 3 ) 
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where f (x) is the function that describes the disturbance. 
Make e << 1, so that the amplitude of the disturbance will 
be small compared with the dimension of the waveduct. 
Assume a solution for this case 

0 = + c 0^ + $2 + . . . + e n $ n + . . . (2.4) 

In the above expression, Oq is as given in equation (2.1) , 
representing the solution for the undisturbed waveduct. 
Since the walls are pressure release, the boundary 

z ~ 0 

z = Z z + e f (x) 



Writing 0 as a Taylor series expansion about the point 
z = l z : 



conditions are: 



0 = 0 at 



0q = 0 at 
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3 4>, 



$U Z + e f (x) ] = $ 0 (£ z ) + e + f(x) ] 



3$ 



1 . 1 ^2 



3 2 4>, 



+ e [ 4-2 + fU) 3 ^ + f f ‘<x) 9 ‘l 

3z j 



+ 



+ £ n [ I 



n ,,3 , . S - 5 $ 

r f J (x) n- j 



j =0 



Jl 



3z 



j X 



+ ... 



( 2 . 6 ) 



With the imposition of boundary conditions, Eq . (2.5) 

then yields the following set of equations: 



3$ 



0 



3z 



f (x) 



3$. 



3z 



f(x) - 



i 3 \ 



21 3z 2 



q ■ f 2 (x) 



n 



? dlxL a3 *n-i 

j=l 2! 3Z 3 



(2.7) 



In the assumed solution, $ and $ n must satisfy the 
homogeneous, linear wave equations 

2 



n 2 A ,_2 13 

□ = (V y — ^ 

C 3t 

2 



) $ = 0 



(2.8) 



n 



,2 1_ d_ 

c 2 3t‘ 



□ ^ = (V ^ -—j) ^ = 0 . 



n 



14 



Substitution of (2.1) into the first equation of (2.7) gives 



= - k cos (k z) cos (cot - k x) 

z z x ' 



. f (x) 



= - k cos(nTr) cos(cot - k x) f(x) . (2.9) 

z x 



Equation (2.9) may be generalized by making it complex. 



±i 



r t \ i / » i(cot-k x) 

= - f(x) k z cos (nTT) e x 



( 2 . 10 ) 



such that 



Re{ $ 



^1 



} = $ 

v 1 



( 2 . 11 ) 



The solution will then become 



= Re{_$ 1 ) 



( 2 . 12 ) 



Define 



F (x) = f (x) e ^x x 



(2.13) 



The insertion of definition (2.13) into Eq. (2.10) produces 



= -e lwt k z cos(nir) F(x) 



(2.14) 
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By use of the definition of the Fourier transform of 



F (x) 



+ CO 



G(s) = / F(x) 



-ixs 



dx 



and its inverse 



(2.15) 



+ c 



F ( x) = / £(s) e dx , 

— OO 



the boundary condition, Eq . (2 . 14 ), becomes : 



±i 



= -e 



icot 



k cos(nTT) ( G(s) 

•z 2 it J — v 



ixs 



ds 



(2.16) 



Use as a trial solution: 



±i 



icot , , . 1 r sm Kz 

■e k cos (nr) — — r T ■ 

z 2tt j sm(K£ ) 

-00 2 



„ , . IXS , 

G(s) e ds . 



(2.17) 



Since the above equation must satisfy Eq . (2.8) [the 

homogeneous, linear wave equation], the quantity K must obey 
the relationship 



K 



[(-) 2 - s 2 ] 1/2 



(2.18) 
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B. SINUSOIDAL UPPER WALL 

Let the upper wall of the waveduct be perturbed in the 
form of a cosinusoide of amplitude where e << 1 . 

(See Figs. 1 and 2) The upper boundary condition is 



$ = 0 at z = & z (l + e cos yx) 



Thus f(x) is specified as 



f(x) - £ z cos(yx) . 



(2.19) 



Inserting Eq . (2.19) into Eq . (2.13) 



F (x) = k cos (yx) e 
1 z 



-ik x x 



iyx , -iyx -ik x 

* I ] • e x 

z 2 



( 2 . 20 ) 



The Fourier transform of (2.20) is: 



k 7 ~ -i(k - Y )x -i(k +y)x • 

G(x) = J [e + e ] e ds 



= ttZ z ISis+k^-y) + 6 (s+k^+y) ] 



( 2 . 21 ) 



and the solution will be: 



l tot . , N 1 r sm Kz 

= ~e k cos (nir) J ■ < rrr 



2? J siiTRiT" \ " £ z [fi (s+k x -Y) + S (s+k x + Y ) 1 



e 1XS ds . (2.22) 
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FIGURE # 1. The Unperturbed Duct 



t+i 




.FIGURE. 2. The Perturbed Duct 
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Using the properties of the 6-function, the integral is 
readily evaluated so that 



. icot , 

$.. = -e k cos nn • 

—1 z z 



. \ 2 1/2 
sm[( c ) (k^ y) z i [cot- (k^+y) x] 

J e X + 

1 • r /^ N 2 ^ , 2 , 1/2 0 

sm [ (— ) - (k x +y) ] Z z 



. r , 0 )* 2 , , . 2 , 1/2 

Sln[( c ) ' (k x" Y) 1 z i[iot-(k -Y)x] 

X 

+ — e 

. r , oo. 2 ,, .2,1/2 0 

sm [ (— ) - (k x -y) 1 



. (2.23) 



Since = Re { } , the first order correction i: 



$, = -k Z cos(nTi) 
—1 z z 



1 J sin K + 

2 \ sin K“ 



cos [cot - (k x +y) x 3 + 
z 



, sin K z r , ,, , , i 

+ sin K-£ cos[(ot ~ (k x~ Y)x U ' 
z 



(2.24) 



where K + = [ (— ) 2 - (k ±y) 2 ] . (2.25) 

V** X 

1 . Validity and Confirmation of Solution 

The above solution is consistent with the results 
presented by N. E. Davis and M. W. Kenyon [1]. The waveguide, 
in that case is perturbed by a traveling surface wave of 
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amplitude A. The upper boundary condition is 4> = 0 at 



7-1 + e F(x,t) , where 

z 

F(x,t) = -l z cos(ftt-yx) and e = 

z 

By making ft = 0 in Eq. (2.12) of [1] and discarding the 
functional dependence in the y-direction, the result agrees 
completely with Eq . (2.22) . 

2 . Properties of the Solution 

In Eq. (2.24) , the propagation constant in the 
z-direction is 



K* = [(f) 2 - (k ±y) 2 ] 1/2 



Since (— ) ^ = k ^ + k ^ , 
c x z 



(2.26) 



+ 



K" 




[1 



k 

(1 ± 2 — ) 
X_ 

(k z / Y ) 2 




(2.27) 



If K 1 



is real and: 



K + Sl = m-fT , m=l,2,..., (2.28) 

the conditions for a resonant situation exist. From 

n 7T 

Eq. (2.2) , k z = — , n=l,2,... . Inserting this into 
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Eq. (2.28) the following expression results, which states 
the relation necessary for a resonance: 



n [ 1 



(1 ± 2 k/y) 

~~ 1 



1/2 



m 



or 



(-) 2 
n 



1 



(1 ± 2 k x /y) 

<Vy > 2 



(2 



The resonance curves are shown in Figure 3 as a plot of 
2 

versus (k /y) . This is in accordance with Eq. (2.16) 

z 

Figure 2 of [1] , and [3] . 



.29) 



k /y 
x' 1 

and 
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III. ANALYSIS 



Using spectral analysis it is possible to find the 
solution to the problem involving any special boundary 
condition related to the upper wall of the waveduct. The 
only thing that has to be known is the Fourier transform of 
the combination 

-ik x 

f(x) • e , (3.1) 

where f(x) represents the function describing the upper wall 
of the waveduct. 

Knowing this, the problem is reduced to solve the 
integral 



T 1 r sm Kz ^ , x ixs , 

I = — rrr— G(s) e ds , 

2 tt J sin K£ — ' 

-oo 2 



(3.2) 



where G(s) is the Fourier transform of (3.1), 



G ( s ) = / f(x) 



-i (k +s) x , 
x dx 



(3.3) 



In the problem treated and solved before, the upper wall was 
described by a very simple function and the evaluation of 
the integral did not pose any special difficulties. But, if 
f(x) is a more complicated function, a careful analysis of 
(3.2) is necessary in order to solve it. 
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A. THE SOLUTION INTEGRAL 

If s is generalized to a complex quantity (s is the 
variable of integration) , a careful study of the behavior of 
the function under the integral sign in (3.2) in the complex 
s-plane is necessary. 

The function to study is: 



f (s) 



. r 2 , 1/2 

sin [ (-) - s ] z 



. r 2 , 1/2 

sin [ (— ) - s ] 



G (s) e 



1XS 



(3.4) 



1 . The Poles 

Assuming, for the moment, that G(s) does not have 
any poles, interest is focused in the ratio of sines in (3.4). 
Since [2] 



sin z = z n ( 1 - 
m=l 



2 2 
m tt 



(3.5) 



f(s) can be written as 



f ( s ) = 



r / Cl) \ 2 2 , 1/2 _ ), 

[ (-) -s ] z n u 

m=l ' 



I(f) 2 -S 2 ] z 2 



2 2 
m -it 



G (s ) e 



1XS 



1 (^) 2 -s 2 ] 1/2 i n i 

C Z , \ 1 «- 

m=l / 



[ (§) 2 -s 2 ] i z 2 



(3.6) 



or simply 
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f (s) 



z 



oo 



(3.7) 



f 2 2 r 2 2, 2 

m -it - [ (-) -s ] z 



'z m=l 



G ( s) e 



1XS 



2 2 r / 2 2. 2 

m 7T - [(-) -s ] z z 



There are no branch cuts, the poles of f(s) are simple and 
located at 



m 



= ± 



- <^> 2 i 



1/2 



(3.8) 



m= 1 , 2 , . . . 



It is possible to arrive at the same result by simply looking 
at the denominator of (3.4). The poles will be located at 



[(|) 2 -S 2 ] 1/2 



Z = ± mu 



(3.9) 



m- 1,2,... . 



Manipulation reveals that Eq . (3.9) is equivalent to 

Eq. (3.8) , but here it is not evident that the poles are 
simple, and that there are no branch cuts. 

From Eqs. (2.2) it is known that 





and 
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n=l , 2 , . . . 



k 

z 



nu 




/ 



Using this in Eq. (3.8) , it can be shown that the location 
of the poles are at 



m 



± [k 



x 



, 7T , 2 , 2 2, , 1/2 

(— ) (n -m ) ] 

z 



n ,m=l , 2 , . . . 



(3.10) 



Looking at this expression, it can bee seen that the poles 
will be located in the real axis if 



, 2 . , tt . 2 , 2 2. 

k x (n “ m } 

z 



(3.11) 



n ,m=l , 2 , . . . 



and on the imaginary axis if the inequality is reversed. 

Let us choose n=l for all that follows. This means 
that only the lowest mode can propagate energy. It can be 
seen that as m increases the associated poles are found 
successively closer to the origin on the real axis until 
above some critical value, m , the poles lie on the imaginary 
axis with distance from the origin increasing with m. The 
first pole, s^, will be located at 

s n = ± k . (3.12) 

lx 
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If it is assumed that m c = 3 , the fourth pole is already 
on the imaginary axis. In Table I and Figure 4 a situation 
like this is shown. 

2 . Effects of Absorption 

If the effects of absorption are introduced, the 
wave equation that governs the propagation inside the duct 
is modified and becomes: 




n 

c 



3 1 



) $ 



n 



= 0 



(3.13) 



n = 0,1,2,... 



where the a n ' s are the coefficients of absorption responsible 
for losses. The a^’s can be assumed as having the same 
magnitude since all the $ ' s have identical form, and the 
propagation constant in the x-direction will become complex. 

It is known that the $'s have the form 

i [ cot-k x] 

sin(k z z)e (3.14) 



Using this in Eq . (3.13) , the value of k is 

x 




ia 



co 

ck 

x 



k 



x 



i& 



where 



3 = 



a 



0) 

ck 

x 



(3.15) 
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TABLE I 



The Poles of f(s) 



n=l , m=l : 



S 1 = f k x 



n= 1 , m=2 



s 2 = ± [k x 2 - 3( x -) 2 ] 1/2 



n=l , m=3 : 



s 3 = * [ k x 2 - 8(^) 2 ] 172 



n=l , m=4 : 



s 4 = ± [k 2 - 15 (-i) 2 ] 1/2 

z 



= ± i [15 (■£”) 2 - k 2 ] 1/2 

z 



n=l , m=5 : 



s 5 = ± i [24 (j —) 2 - k x 2 ] 1/2 

z 



n= 1 , m= 6 , 7 , 



• r/2,,,77.2 , 2, 1/2 

s m = ± 1 [(m -1) (— ) - k x ] 

z 
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a 

t 



if 



A s <> 



A 4 <j> 




IU 



<>/i, 



ra c =3 




FIGURE 4. The Poles of f( s ) 
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is the total absorption coefficient in the x-direction and 
will be considered very small. 

It is necessary to know what is the effect of the 
introduction of absorption in the location of the poles of 
f(s), especially the ones located in the real axis. If s-^ 
is located at s, = ±k , it can be seen that the poles 
in the negative part of the real axis will be moved up and 
the poles located in the positive part will move down 
(Figure 5) . We may now choose a contour in order to solve 
the integral. 

B. THE CONTOUR INTEGRAL 

The integral to be solved is 

CO 

I = J f (s) ds (3.16) 

— co 

where f(s) is the expression in (3.4). 

1 . Choice of the Contour 

The integral in (3.16) is very similar to the one 
Isakovich used in his work [4]. It reduces to 27 Ti times 
the summation of the residues at the poles of f (s) . 

The contour to consider will consist of the segment 
[-R,+R] on the real axis and the upper semicircle C R , as 
shown in Figure 5. 

Therefore 

°° +R 

I = / f (s) ds = lim j f(s) ds = 2ni £ Res^ f( s ) • (3.17) 

-co R-><» -R m 
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I 

! 




FIGURE 5. 



Effects of Ahsorp 
°f the Poles of f 



on the Location 
Contour . 
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Simplify Eq. (3.4) by defining H(s) as: 



H(s) 



. r , to, 2 2 ■, 1/2 

sin [ (— ) - s ] z 



• r / (jO\ 2 2 ■* f. 

sin [ (— ) - s ] £ z 



(3.18) 



For the case where G(s) does not have any poles the 
integral will be: 

co ixs 

I = / f (s) ds = 27ri £ G (s^) e m Res{H(s m ) } (3.19) 

-CO m 



On the other hand, if the Fourier transform of F(x) has poles 
of its own. at s = s , where 

g 



gWT * 0 ' (3 - 20) 

- g 

the expression to evaluate the integral is modified as 
follows : 



00 ixs 

I = / f(s) ds = 2 Tri £ G(s ) e m Res{H(s ni )} 

_co m 

ixs 

+ 27ri 'l H (s ) e g Res{G(s )} (3.21) 

g 9 9 



In either case, the poles of interest will be the ones 
included within the contour chosen. 
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2 . 



Definition of the Transform of F(x) 



In order always to be able to close the contour in 
the upper half plane, the definition of the Fourier transform 
of F(x) must be altered. 

We define 

CO # 

G ( s) = / F(x) e T1SX dx 

— 00 



with inverse 



F(x) = ~ / G(s) e ^ XS ds 

— 2 7T ; — 

— CO 

The upper sign of the exponent will represent the 
positive x-part of the waveduct with waves standing or 
propagating to the right and the lower sign will represent 
the negative x-part of the waveduct with waves standing or 
propagating to the left. 
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IV. SPECIAL CASES 



In order to simplify the problem and apply this theory 
to some specific cases we shall impose some conditions on 
the problem. Since 

(^) 2 = k 2 + k 2 , (4.1) 

C X z ' 

in the duct, the cut-off frequency will be 



u = c k (4.2) 

The sound source that is going to excite the duct will have 
a monofrequency of 

a) = A c k , (4.3) 

z 

where 1 < A < 2 (4.4) 



to assure propagation only in the lowest modes. Under this 
condition and knowing that 



k 



z 




(4.5) 



(only lowest mode allowed in z, n=l) the location of the 
poles of f(s) will be at 
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( 4 . 5 ) 



It can be seen that only one pole will be on the real 
axis (s-^) , and that all others are already on the imaginary 
axis . 

The general expression for the location of the 
poles of f(s), assuming G(s) has no poles is: 



s 



m 



± [k 2 + (-— -) 2 (1 - m 2 ) ] 1/2 

x Z 



( 4 . 6 ) 



and the rn residue is 



Res 



m 




z 



2 



( 4 . 7 ) 



evaluated at s = s 



m 



This will become, after manipulation 




z 




( 4 . 8 ) 



m 



1,2 



f • • • 
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th 

The effect on acoustic velocity potential due to the in 
pole is [using Eq. (2.17)] 

el “ t k z TZ ' 2l,i l Res ra ’ (4 ' 9) 

If G(s) has poles of its own, (4.9) will be modified 
and the effect on the velocity acoustic potential will be 



I 

m 

+ 



A. THE 6-FUNCTION DISTURBANCE 

Define f (x) , the function describing the wall disturbance 
as 



iwt , 
e k 



_ 1 _ 

2 TT 



2 ni 



Res 



m 



±ixs 

l H ( s ) e g Res_-j (G(s_) } 

g g 



(4.10) 



f (x) = -£ z 6 ( x) 



(4.11) 



Then 



-ik x 

F (x) = -£ 6 (x) e 

— z 



(4.12) 



and 



G (s) = -£ z 



(4.13) 



with no poles . 
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The poles of f (s) will be located according to Eq. (4.6) . 
The first pole, the only one in the real axis will be at 



s, = -k (4.14) 

lx 

and the second one, on the imaginary axis will be at 

s 2 = + 1 [ ( £-) 2 ( 2 2 - 1 ) - k 2 ] 1/2 (4.15) 

z 

= + iP 



where P is defined as 



P = ~ k x 2 ] 1/2 

z 



(4.16) 



The overall situation is shown in Figures 6 and 7. The 
residue of f(s) at s^, Res-^(s^) will be, using Eq. (4.8), 



k z ?ik x x 

Res^(s^) = — sin k^z e 

x 



(4.17) 



The acoustic potential due to this residue will be calculated 
using (4.9) , and 



— i{R es i (S]_) I = 



k i(o)t+k x) 

z , . x 

— sin k z i e 
k z 

x 



(4.18) 
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FIGURE 6 



The 6-Function Disturbance 




FIGURE 7. The Poles of Interest 



(4.19) 



Since 



$ 



1 



Re ( jL x ) 

2 



x 



, in this case 



sin k z sin(ojt+k x) 
z x 



The above equation represents two scattered traveling waves 
emanating from x=0, one going in the positive x-direction and 
the other in the negative x-direction. 

The residue of f(s) at S 2 , Res^{s n ) , will be, using 
Eq. (4.8), 



Res2 ( S 2 ) 



2k 

z 

-1 — sm 



2k z 
z 



+ xP 
e 



(4.20) 



with a 



corresponding acoustic velocity potential of 

2k ^ 

[Res 2 (S 2 ) ] = — jr~ sin 2k z z e^ wt e + xP 



(4.21) 



This represents two disturbances in the acoustic field, 

decaying in the ± x-direction . It can be seen that far away 

from the wall disturbance their effect will be negligible. 

Very close to the wall disturbance, x=0, the effects will 

be considerable and depend on the size of P m , diminishing 

as P_ increases with m. P_ is defined as 
m m 

•n _ r , 2 . w 7i N 2 , 2 , 1/2 _ 0 „ 

Pj^ [ (m — 1) ( ^ ) ~ ^ x ' ^ ^ f * * * * 

z 

However, even at great distances from the irregularities 
in the wall the acoustic field is affected by [Res-^ ( s-^) ] . 
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The total field, discarding the effects of decaying 



disturbances as in (4.21) is: 




k 



2 



= sin k z[cos(wt-k x) 
z x 



z 



sin (ojt+k x x) ] 



(4.22) 



k 



x 



Note that the traveling perturbation is 90° out of phase with 
the initial wave. 

B. THE DECAYING WALL 

Define f(x), the function describing the wall disturbance 
as 



f(x) = -Z z e 



ax 



(4.23) 



Then 



F(x) = -Z z e 



ax 




(4.24) 



e 



and 



G ( s ) 




(4.25) 



I a I [1 + (- + k ) 2 ] 

II a x 



with poles located at 




(4.26) 
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The poles of interest of f(s) will be located at 



s 2 = +iP ( J .27) 

s = a (-k + i) 

q x 



The overall situation is shown in Figures 8 and 9. Since 
G(s) has poles, Eq . (3.21) must be used to evaluate the 

integral . 

The residue of f(s) at s^, Res^(s^) will be 



2 k z 

Res, (s,) = — 7 — sin k z 
11 a k z 

x 



ik x 
x 



2 , . 2 , ..2 

a + k^ (a-1) 



(4.28) 



with a corresponding acoustic velocity potential of: 



[Res-^ (s 1 ) ] 



i 2 



k 2 , , 2 , _ . 

x a + k (a- 1 ) 
x 



tt sin k z 
2 z 



i ( tot +k x x) 
e 

(4.29) 



representing two traveling waves emanating from x= 0 , going 
in opposite directions, and 



2k 



= Re { } = 



: 2 , , 2 , , . 
x a + k (a- 1 ) 
x 



~rr sin k z sin(a)t+k x) 

^ Z X 

(4.30) 
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FIGURE 8. The Decaying Wall 



1 




FIGURE 9. The Poles of Interest 
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The residue of f(s) at S 2 > Res 2 (s 2 ) will be 



4k. 



ReS2(s2) = 



.iP 



iPa[l + (— + k x ) ] 



7T— sin 2k z e 
2 , z 



+Px 



(4.31) 



with a corresponding acoustic velocity potential of 



4k 



ilt Res 2 (s 2 )] = ~ AS 



Pa[l + (^-) + k )‘] 

Cl X 



lU)t +Px 

rr~ sin 2k z e e 
2 , z 



(4.32) 



This represents two standing disturbances decaying with 
x and negligible at great distances from the origin. 

The residue of f(s) at s will be computed using (3.21) 
and is 



± 1XS 

H (s ) e g Res{G(s ) } 

y y 



(4.33) 



or 

k x 

sin /B+iC z +ax + x . „ 

e e 1 i 

sin /B+iC 

where : 

2 2 2 2 
B = k^+k^(l-a) + eT 
z x 

C = 2a 2 k 

x 



(4.34) 



(4.35) 
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The corresponding acoustic velocity potential is: 



[Res (s ) 3 
1 g g 



sin /B+iC z 
Z Z sin /B+iC 1 

z 



+ax 

e 



i (cut+k x x) 
e 



(4.36) 



This corresponds to two traveling-decaying waves (both 
in x and z directions) going in opposite directions from 
x=0. At large distances from the origin its effect will be 
practically nil. 

Combining results, only valid at x >> 0 or x << 0 , 
we have: 



$ = + e 



= sin k z [cos(cot-k x)' 

Z X ' 



2k ^ 

z a 

k 2 , 2 , , N 

x a + k ( a-1) 
x 



-j sin ( oot+k x ) ] 



(4.37) 



It can be noticed, again, that the traveling disturbance 
is 90° out of phase with the original propagating wave. 
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V. RESULTS AND CONCLUSIONS 



The results obtained show that the effects of very small 
irregularities in the physical boundaries of the duct can 
have appreciable importance in the waveduct sound propagation 
even at very large distances from the disturbances. In the 
case of the sinusoidal wall, resonances can exist, which 
could have the same order of magnitude of the original 
propagating wave. These resonances do not appear in the 
two other special cases studied, the 6-function disturbance 
and the decaying wall. It seems that the theory works and 
can be applied to other types of corrugations , once the 
Fourier transformation of the function describing the wall 
disturbance is known. However, complications can arise 
depending on the complexity of the above mentioned Fourier 
transform, making the solution integral difficult to evaluate. 
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